Finite bias Cooper pair splitting 
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In a device with a superconductor coupled to two parallel quantum dots (QDs) the electrical 
tunability of the QD levels can be used to exploit non-classical current correlations due to the 
splitting of Cooper pairs. We experimentally investigate the effect of a finite potential difference 
across one quantum dot on the conductance through the other completely grounded QD in a Cooper 
pair splitter fabricated on an InAs nanowire. We demonstrate that the electrical transport through 
the device can be tuned by electrical means to be dominated either by Cooper pair splitting (CPS), 
or by elastic co-tunneling (EC). The basic experimental findings can be understood by considering 
the energy dependent density of states in a QD. The reported experiments add bias-dependent 
spectroscopy to the investigative tools necessary to develop CPS-based sources of entangled electrons 
in solid-state devices. 
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The electrons of a Cooper pair in a conventional super- 
conductor form a spin singlet, which might be exploited 
as a naturally occurring on-chip source of spin-entangled 
Einstcin-Podolsky-Rosen (EPR) P, i] electron pairs, if 
the electrons can be separated coherently. This Cooper 
pair splitting (CPS) can be understood as inverse crossed 
Andrccv reflection and was initially searched for in metal- 
lic nanostructures with tunnel contacts 0-0] • However, 
other processes like elastic co-tunneling (EC) make the 
detection of CPS difficult. Electron-electron interactions 
are relevant in these processes, but can not be tuned by 
external means in these structures. It has been suggested 
to use electrically tunable quantum dots (QDs) coupled 
to a superconducting lead to obtain such tunability [1, Q . 
Recently, first transport characteristics at zero bias were 
reported for InAs nanowire (NW) and carbon nanotubc 
QD devices, which are explained by CPS 

Bias-tunable non-local resistance in a metallic struc- 
ture was attributed to the excitation of different modes 
of the electromagnetic environment by CPS and EC 
0, This mechanism, however, is difficult to con- 

trol experimentally. Here we report finite bias differen- 
tial conductance measurements on InAs nanowire devices 
similar to those in Ref. [l^l. We show that the conduc- 
tance through one QD is not entirely due to local pro- 
cesses, but also due to non-local higher order tunneling, 
consistent with a simple picture of CPS and EC. Our ex- 
periments establish an additional parameter and a proce- 
dure to characterize a generic Cooper pair splitter. Since 
the relevant processes depend differently on the bias, the 
relative rates can be tuned by external means. In particu- 
lar, we show that the energy dependent effective density 
of states (DoS) due to the QD levels is crucial for the 
reported effects. 

A colored SEM image of a generic CPS device is shown 



in Fig. 1(a). An InAs wire [ij] is contacted by a super- 
conducting strip in the center (S) and two Ti/Au normal 
metal leads (Nl, N2), which define two quantum dots 
QDl and QD2. The two top-gates (gl, g2) are strongly 
decoupled from each other, so that the QD levels can be 
tuned individually, while the highly-doped Si wafer sepa- 
rated from the device by 400 nm thermal oxide serves as a 
global back-gate. More details on the sample fabrication 
can be found in [l3| . All experiments were performed at 
the base-tcmpcrature of T « 20 mK. 

Also shown in Fig. 1(a) is the measurement schematic: 
an ac voltage of U^c = 10 /iV at 77 Hz is applied to 
the superconductor at zero dc potential. A dc voltage 
C/n2 is applied to contact N2 of QD2, while a home-built 
current- voltage converter (gain 10^ V/ A) and a standard 
lock-in amplifier are used to measure the ac current Ii 
in contact Nl from QDl. A 10 ^F capacitor C is used 
to decouple possible IV-converter off-set voltages from 
the device, while keeping the dc potential at zero. Fig- 



(b) 




Nt 



N2 



QD1 QD2 
DPT 



FIG. 1: (Color online) (a) Colored SEM image of a typical 
device and measurement circuitry, (b) Schematics of direct 
Cooper pair tunneling (DPT) with an ac voltage dU imposed 
on QDl and QD2 and a finite bias (7n2 across QD2. 
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FIG. 2: (Color online) (a) Characterization of QDl: differen- 
tial conductance Gi as a function of Uni and f/gi at B = 0. 
The cross section shows two small features related to the su- 
perconductor gap A. (b) Gi as a function of Un2 and Ugi at 
5 = 0. (c) AGi = Gi(f/gi,C/N2) - Gi(f7gi,C/N2 = -ImV) 
derived from the data in (b) for the same parameter range, 
(d) AGi for the same experiment at B = 200 mT. 



urc 1(b) shows the energy diagram of such a device. 

To characterize QDl, the differential conductance 
Gi = dli/dU is plotted in Fig. 2(a) as a function of 
the gate voltage [/gi and the bias applied to Nl, J/ni. 
We observe well-defined Coulomb blockade oscillations 
consistent with the formation of a QD IJ]. We label 
the resonances in the figure as R1-R3 and the electronic 
states in-between as Sa and Sb, from left to right. From 



this plot we extract various sample parameters, e.g. the 
superconducting gap A « 130 /icV (cross section at 
C7gi = —0.146 V), or an addition energy of ^ 3meV and 
the tunnel coupling of F »i 500 ^.eV for Sb. In addition, 
a Kondo resonance [l5| in Sb can be observed at zero 
bias in the normal state data (not shown). A similar 
characterization of QD2 in the parameter range of the 
experiments (not shown) exhibits only small variations 
in the conductance. In the discussion below we therefore 
assume a resonance in the effective DoS of QDl and a 
constant DoS for QD2, as illustrated in Fig. 1(b). 

The dominant component of Ii is due to local trans- 
port processes between S and A^i. As an example, direct 
Cooper pair tunneling (DPT), also known as Andreev re- 
flection, is illustrated in Fig. 1(b): applying the voltage 
dU to S as in the experiment is equivalent to applying 
—dU to Nl and N2. An electron at energy E and one 
at —E are injected from Nl to form a Cooper pair in 
S. Because S and A^i are at the same dc potential, all 



local transport processes are independent of the voltage 
U^2 on N2, if charge imbalance jl6[ and local heating 
can be neglected. Therefore only 'non-local' transport 
processes, which comprise correlated coherent tunneling 
of electrons through QDl and QD2, can lead to a depen- 
dence of Ii on ?7n2- 

Figure 2(b) shows Gi as a function of the top gate 
voltage ?7gi of QDl and the bias ?7n2 applied to QD2. 
On this scale the conductance seems independent of U^2 ■ 
We assume that all non-local processes become ineffec- 
tive at [/n2 ^ A/e and therefore plot in Fig. 2(c) the 
deviation of Gi from a high-bias value, AGi(C/gi, [/n2) — 
Gi([/gi, [/N2)-Gi(t/gi, J7n2 = -1 mV). We note that the 
colorscale is adjusted to white representing AGi ~ 0. For 
the gate voltage Ugi tuned to resonance Ri, AGi has a 
positive maximum at zero bias with AGi/Gi ~ 15%. 
The extent of this maximum is roughly eJ7N2 ~ A and 
eaUgi ~ r {a: lever arm of top gate gl). With i7gi 
slightly off-peak, a minimum in AGi occurs on both sides 
of resonance Ri, but not centered around J7n2 = 0. These 
features will be examined in more detail below. Similar 
features occur on the resonances 2 and 3, however with- 
out minima where state B is involved. 

In a control experiment the superconductivity is sup- 
pressed by an external magnetic field of 200 mT paral- 
lel to the Al strip. The corresponding plot of AGi is 
shown in Fig. 2(d). No clear structure can be discrim- 
inated, though the scale is considerably smaller than in 
Fig. 2(c). We therefore conclude that the features in 
Fig. 2(c) are due to the superconductor and non-local 
transport processes, which become unlikely if S is in the 
normal-conducting state. 

Figure 3 shows cross sections AGi([/n2) of Fig. 2(c). 
The corresponding top-gate voltages Ugi are near reso- 
nance 1 and indicated on the zero bias top-gate sweep 
Gi(C/gi) in Fig. 3(a). For gate voltages far from the res- 
onance the bias dependence of AGi shows a small but 
pronounced minimum at zero bias, as shown in Figs. 3(b) 
and (g), while on resonance a strong peak is observed at 
zero bias, see Fig. 3(d). Zero bias across QD2 is indi- 
cated by dashed vertical lines in all plots and is slightly 
offset due to a small input offset-voltage from an IV- 
converter mounted on lead N2. Slightly off-resonance 
the bias dependence is asymmetric with respect to zero 
bias: in Fig. 3(c) AGi exhibits a minimum at negative 
and a maximum at positive bias, while on the other side 
of the resonance a maximum can be found at negative 
and a minimum at positive bias, see Figs. 3(e) and (f). 
We note that the maxima and minima near the resonance 
are not at zero bias, whereas far off the resonance only a 
minimum at ?7n2 = is found. 

The non-local signals decay with increasing tempera- 
ture, as shown in the inset of Fig. 3(a) for the conduc- 
tance minimum at Ugi = —0.1485 V, essentially the gate 
position of Fig. 3(g). The amplitude of AGi decreases 
monotonically and disappears at T sa 175 mK. Up to this 
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FIG. 3: (Color online) Bias dependence of AGi for a series 
of top gate voltages Ugi. The latter are indicated in (a). 
The dashed curves are derived from the model discussed in 
the text. Vertical lines indicate I7n2 = 0, including a small 
offset from the IV-converter on N2. The inset in (a) shows 
the temperature dependence of the the minimum in (g) with 
a black line as a guide to the eye. 



temperature we have found no significant change in the 
superconductor gap. We therefore conclude that A is not 
the limiting energy scale in this problem, reminiscent of 
the CPS observed in [l^l at zero-bias. 

A qualitative understanding of the experiments can be 
gained by expressing the CPS and EC tunneling rates by 
energy dependent effective density of states Di{E) and 
D2{E) in QDl and QD2, incorporating the respective QD 
transmissions. The processes are illustrated in Fig. 4, for 
which we assume a resonance feature for Di and a con- 
stant for QD2, i.e. D2{E) = £>2- We neglect changes in 
the DoS by virtual tunneling processes 18| and electron- 
electron interactions in the QDs and the leads. 

The probability that an electron from QDl at energy 
E and one from QD2 at energy ^E form a Cooper pair 
in S (inverse CPS) depends on the DoS as 



K McPs{E)Di{E)D2{~E)f{E,ti^)f{~E,ii2), (1) 

with the respective electrochemical potentials /ii_2 in the 
leads, measured relative to the = in the super- 
conductor. Mcps represents the transition probability. 
For CPS the states of the QDs have to be empty ini- 
tially and the Fermi functions /(±_E, are replaced 
by 1 — f{±E, ^1^2)- We note that generally we have in- 
verse CPS for C/n2 < and CPS for C/n2 > 0. The situa- 
tion for EC is shown schematically in Fig. 4(b). Since the 
process is elastic, its probability at the energy E scales 
as 

cx MEciE)Di{E)D2{E) [/(£;, ^^2) - f{E,^ll)] . (2) 



Summation over all energies for periodically varying [/ni 
and U^2 results in the non-local (nl) differential con- 
ductance, AGi = dli^ni/dU. Assuming constant Mcps, 
A/ec and D2, as well as T — > and edU << F, the fol- 
lowing intuitive expressions for the contributions of CPS 
(Gcps) and EC (Gec) to AGi = Gcps + Geo can be 
found: 



Gcps - -j^McpsD2[Di{0) + D^i+eUm)] (3) 



G 



EC 



-A/ec^2 [Di{0) - Di{-eUf,2)] 



(4) 



Both processes have a component from the DoS at E' = 0. 
In addition, CPS has a positive contribution from Di at 
+eU , while EC has a negative contribution from — eC/. 
This can also be understood from the diagrams in Fig. 4. 
We note that in our measurement scheme the differential 
conductance due to EC would be zero for energy indepen- 
dent transmissions, e.g. in metallic structures, and the 
contribution of CPS would be constant and independent 
of the bias applied to N2. 

First we discuss the case when QDl is tuned to a sym- 
metric resonance, i.e. Di{E) — Di{~E) and l?i(0) > 
Di{E ^ 0), as depicted in Fig. 4(a) and (b). From 
Eqs. ([3]) and (jlj one finds that Gcps has a maximum 
and Gec a minimum with Gec = at [/n2 = 0. Both 
signals are positive on a resonance. If the two processes 
had the same probability, Mcps = M-^c, the variations 
of Gcps and Gec would cancel exactly and we would 
not expect any changes in AGi as a function of Un2- We 
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FIG. 4: (Color online) Schematics of CPS and EC at finite 
bias on QD2 and an ac voltage dU imposed on both QDs. (a) 
and (b) show the situation when QDl is tuned to a resonance, 
while (c) and (d) depict an off-resonance case. 
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therefore conclude that CPS and EC can obtain different 
relative weights and the comparison with the experiment 
in Fig. 3(d) suggests that on a resonance CPS is the dom- 
inant process, i.e. M^ps ^ -^eCj independent of 

An off-resonance situation is illustrated in Figs. 4(c) 
and (d). If the resonance is shifted in energy to E' = i?o = 
—eaAUgi by top-gate 1, Eqs. ([3]) and (|4]) predict two fea- 
tures: 1) Since Di{Ea) > Di{0) > Di{-Eo), one finds 
for eJ7N2 = Eq that Gcps has a positive maximum and 
dominates AGi. 2) Similarly, one finds for eU^2 = —Eq 
a negative minimum in both, Gec and AGi, the latter 
because Gcps is negligible. Such features are observed in 
Fig. 3(c) and (e), where a minimum and a maximum oc- 
cur as a function of J7n2, with reversed order depending 
on the gate voltage relative to the resonance. 

The results of numerical evaluations of Eqs. (2) and 
(3) with a Gaussian for Di and the adjustable parame- 
ters Mcps, Mec and an offset, leads to the dashed curves 
in Figs. 3(b)-(g). The data are reproduced qualitatively 
near the resonance and we find that CPS is favored on 
resonance, while slightly off-resonance EC obtains a sim- 
ilar strength as CPS. These findings might be related to 
a reduction of CPS for energies close to A, while EC is 
increased due to the energy-dependence of the second or- 
der tunnel matrix elements. Both rates diminish above 
A due to the creation of excitations in S. 

This model is very rudimentary and neglects electron- 
electron interactions in the leads and on the QDs, which 
might be crucial to understand details in our data. A 
strong indication for this is that in all Coulomb blockade 
regions a minimum occurs at zero bias, see Figs. 3(b) 
and (g), which can not be reproduced in the model. The 
dashed curves in these figures are for Mec ^ Mcps and 
show a minimum adjusted to the dip in the data, but at a 
wrong bias. A dip at zero bias, however, is in qualitative 
agreement with a zero-bias anomaly due to dynamical 
Coulomb blockade at low transmissions @. In addition, 
in state Sb, for which we observe Kondo correlations in 
the Coulomb diamonds, the minima in AGi expected 
near the resonances are missing or very weak, while they 
are clearly observed for the other states. This might be 
due to a competition between the superconducting and 
the Kondo correlations, which is not accounted for in our 
model. 

In summary, we have reported finite bias measure- 
ments on an InAs nanowirc quantum dot Cooper pair 
splitter. Our results show that the energy dependent 
transmission due to the QDs has strong effects on non- 
local processes in such systems. We show that finite- 
bias spectroscopy is useful to identify non-local processes 
and find that Cooper pair splitting can be the dominant 
process on a QD resonance. Off-resonance we find a pat- 
tern consistent with EC being of similar strength as CPS. 
These findings are interpreted in a simple model based 
on energy dependent effective density of states. However, 



this model does not account for the relative strength of 
the processes, nor for features in the Coulomb block- 
ade regions or where Kondo correlations are relevant. 
We tentatively attribute the latter findings to electron- 
electron interactions, which should be studied in greater 
detail to understand the involved mechanisms and find 
means to exploit them in an on-chip source of entangled 
electrons. 
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